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It is well known that the result of any phase measurement on an optical mode made using linear 
optics has an introduced uncertainty in addition to the intrinsic quantum phase uncertainty of the 
state of the mode. The best previously published technique [H. M. Wiseman and R. B. Killip, Phys. 
Rev. A 57, 2169 (1998)] is an adaptive technique that introduces a phase variance that scales as 
n~^"^ , where n is the mean photon number of the state. This is far above the minimum intrinsic 
quantum phase variance of the state, which scales as n~^. It has been shown that a lower limit to 
the phase variance that is introduced scales as ln(n)/n^. Here we introduce an adaptive technique 
that attains this theoretical lower limit. 

PACS numbers: 42.50.Dv, 03.67.Hk, 42.50.Lc 



I. INTRODUCTION 

The phase of an electromagnetic field cannot be mea- 
sured directly using linear optics and photodetectors. 
Rather than measuring phase directly, phase measure- 
ment schemes rely on measuring quadratures of the field 
and inferring the phase from these measurements. In 
a typical experimental implementation, the mode to be 
measured is passed through a 50/50 beam splitter, in or- 
der to combine it with a much stronger local oscillator 
field. The difference photocurrent from the two output 
ports of the beam splitter yields a measurement of a par- 
ticular quadrature. 

The standard technique for measuring a completely un- 
known phase is heterodyne detection, where all quadra- 
tures are sampled with equal probability. This is achieved 
by using a local oscillator field with a frequency slightly 
different from the signal's frequency, so its phase changes 
linearly with respect to the phase of the signal. More ac- 
curate phase measurements can be made using the ho- 
modyne technique, where the local oscillator phase is 
$ = (/3 -f 7r/2, with if the phase of the signal. The prob- 
lem with this is that it requires initial knowledge of the 
phase of the signal, and so is not a phase measurement 
in the strict sense. 

To maintain the unbiased nature of heterodyne phase 
measurements but obtain the increased sensitivity of ho- 
modyne measurements, an adaptive dyne technique can 
be used ||, ||, Here "dyne" detection is used to 
mean photodetection using a strong local oscillator at 
a beam splitter. The idea behind adaptive phase mea- 
surement schemes is to use the information gained so far 
during the measurement to estimate the system phase ip. 
This is then used to adjust the local oscillator phase $ 
to approximate a homodyne measurement as above. 

The apparatus for performing these measurements is 
shown schematically in Fig. |^. The signal and a local 
oscillator with amplitude (3 are combined at the beam 
splitter and the outputs are measured with photodetec- 
tors. The outputs from the photodetectors, and 
(57V_, are subtracted and then fed into a digital signal 



processor that uses these measurements to estimate the 
phase of the system, and adjusts the phase of the local 
oscillator via an electro-optic phase modulator. The sig- 
nal is shown here as from a cavity with a half-silvered 
mirror, as this is what is considered in the theory in Sec. 

iv. 

The signal of interest is the difference between the pho- 
tocurrents at the two ports. We therefore define the sig- 
nal as 



I{t)5t 



5N+ - SN^ 
/3e-*/2 



(1.1) 



Here we have used units of time such that the decay 
constant of the cavity is unity. We divide by a factor of 
e~*/^ because this is the square root of the mode function 
for the signal. We can take account of signals with more 
general mode functions u{t) in a similar way ||^. 

When we take the limit of very large local oscillator 
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FIG. 1: Diagram of the apparatus for making an adaptive 
phase measurement. The signal from the cavity is combined 
with the local oscillator field at a 50/50 beam splitter (BS) 
and the outputs are detected by photodetectors (PD). The 
signals from these photodetectors are subtracted, and the 
difference signal is processed by the digital signal processor, 
which determines a phase estimate and adjusts the electro- 
optic phase modulator (EDM) accordingly. 
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amplitude and small time intervals St, we find that 

I(v)dv = 2Re(a^e-'*(^))dw + dW(v), (1.2) 

where v is time scaled to the unit interval and is 
the scaled mean amplitude of the system (for which the 
S superscript stands). The systematic variation in the 
coherent amplitude with time due to the mode shape is 
scaled out. This scaling is explained in more detail in 
Sec. The final term dW{v) is an infinitesimal Wiener 
increment such that {dW{v)'^) — dv 

It can be shown |^ that just two complex numbers 
are necessary to encapsulate all of the relevant informa- 
tion in the photocurrent record up to a given time. These 
are 



A, 



I{u)e"^^'''>du, 







Jo 



(1.3) 
(1.4) 



For convenience, we often replace By by a third complex 
number defined in terms of Ay and By , 



Cy A,. 



ByAl 



(1.5) 



Generally the best estimate of the phase at time v is 
arg(Ci,) [||. The subscripts are omitted for the final val- 
ues {v — 1). 

In adaptive measurement schemes the phase of the lo- 
cal oscillator is generally taken to be 



(1.6) 



where <^(w) is the estimated phase of the system at time 
V using the measurement results Ay and Cy. There are 
a number of possible phase estimates, giving different 
adaptive schemes. For the mark I scheme ||, both the 
running phase estimate (p{v) and the final phase estimate 
are taken to be &T:g{Ay). This is better than heterodyne 
measurements only if the field is very weak j|l|, ^. For 
the mark II adaptive phase measurements |^ the best 
phase estimate arg(C) is used at the end of the measure- 
ment, but for the intermediate phase estimate avg{Ay) is 
used. This is better than heterodyne measurements for 
all field strengths. If arg(C^) is generally the best phase 
estimate, it is apparent from Eq. ( |l.5| ) that arg(A^) will 
only be the best phase estimate if By is negligible (as it 
is in the case of heterodyne measurements) . For adaptive 
phase measurements By does not vanish and axg^Ay) is 
generally a much worse phase estimate than aig{Cy). 

This raises the question of why this relatively poor in- 
termediate phase estimate is used. There are two main 
reasons for this: (i) it is possible to obtain direct analytic 
results for this case, whereas using a better intermediate 
phase estimate requires numerical evaluation; (ii) the ap- 
paratus required to implement this method is much sim- 
pler than that required for a better intermediate phase 
estimate. 



Even with the relatively poor intermediate phase es- 
timate, the mark II adaptive scheme introduces a phase 
variance of just \{n^)~^'^ , a good improvement over the 
heterodyne result of ^{n^)~^ . Here is the mean pho- 
ton number of the field being measured, and the actual 
measured phase variance is the introduced phase variance 
plus the intrinsic phase variance. The intrinsic phase 
variance for a state of mean photon number can be as 
small as of order (n"^)^^ |^. This is far smaller than 
the introduced phase variance, so the latter is what lim- 
its the accuracy of phase measurements. Although the 
mark II results are far superior to the standard result of 
heterodyne detection, it is still possible to improve on the 
mark II result, and it is shown in Ref. Q that a theoreti- 
cal lower limit to the phase variance that is introduced by 
an arbitrary phase measurement scheme (based on linear 
optics and photodetection) is j\n{n^) x (n"^)"^. 

In improving on the mark II result, the obvious thing to 
do is to use a better intermediate phase estimate. It turns 
out that using the best phase estimate arg(Ci,) actually 
gives a worse result than the mark II case, for reasons 
that we will explain later. The phase estimates that we 
consider in this paper are therefore intermediate between 
arg(yl,j) and the best phase estimate: 



0{v) = arg(C^'(-)A:(")). 



(1.7) 



It is possible to obtain a marked improvement over the 
mark II case by using constant values of e. We show 
in Sec. ^ that a scaling of roughly {n^)~^-^^ is possible. 
One drawback is that the value of e required depends on 
the photon number. 

We can obtain an even better result if we allow e to 
have a variation in time, and we show in Sec. ^ that we 
can obtain phase estimates very close to the theoretical 
limit if we use 



e{v) 



\By 



Cv 



1 



(1.8) 



This expression does not explicitly depend on the photon 
number. This method works best if the phase estimates 
are updated in discrete time steps, and the magnitude 
of the steps depends weakly on the photon number. A 
more serious problem with this method is that it tends to 
produce values of that are too close to 1. This means 
that final phase estimates with an error close to tt occur 
sufficiently frequently to make a significant contribution 
to the phase uncertainty. We will show how this problem 
can be corrected. 

The paper is structured as follows. In Sec. || we red- 
erive the ultimate theoretical limit to phase measure- 
ments of Ref. Q . This is necessary to unde rstand how the 
improved feedback algorithm of Eq. (1.7) can approach 
the theoretical limit, as explained in Sec. III. In Sec. IV 



we derive the results necessary for a numerical simula- 
tion of this algorithm, and in Sec. ^ present the results 
of those simulations. The problem of infrequent results 
with large errors is identified in Sec. ^ and a solution 
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proposed and evaluated in Sec. |VII[ W e conclude with a 
summary and discussion in Sec. |VII]|. 



II. THE THEORETICAL LIMIT 

In order to understand how to attain the theoretical 
limit, we must first understand the reason for the theo- 
retical limit. It can be shown |6| that the probability of 
obtaining the results A, B from an arbitrary (adaptive 
or nonadaptive) measurement is 

P{A,B)(fA(fB = Ti[pG{A,B)](fA(fB, (2.1) 

where p is the state of the mode being measured. Here 
G{A, B) is the POM (probability operator measure) for 
the measurement, and is given by 



G{A,B)=Q{A,B)\i;{A,B)){i:{A,B)l 



(2.2) 



where Q{A,B) is what the probability distribution 
P{A, B) would be if p were the vacuum state |0)(0|, and 
|V'(^, B)) is an unnormalized ket defined by 



|V(A,B)) =exp 



2 ^ ' 



|0). 



(2.3) 



This is proportional to a squeezed state 



exp 



'-B{a)f-Aa^ 



|0) = (l-|S|2)-'/%xp(Aa72)|a,0, 

(2.4) 



2 

where 

|q!, ^) — exp(Q!a^ — a*a) exp 



and the squeezing parameters are 
A + BA* 



1- |S|2 ' 
Batanh|S| 



B 



|0), 
(2.5) 

(2.6) 
(2.7) 



where atanh is the inverse hyperbolic tan function. In 
terms of these the POM is given by 



G{A,B)^Q'{A,B)\a,0{a,£,l 



(2.8) 



where 



Q'{A, B) = Q{A, B) (1 - |B|2) exp [Re {Aa*)] . 

If the system state is pure, p — \'ip){tp\ and the probability 
distribution is given by 



P{A,B) = Q'{A,B)ma,i)[^ 



(2.10) 



For an unbiased measurement scheme the probability 
distribution for the phase resulting from this equation 



depends entirely on the inner product between the two 
states, and not on Q'{A,B). To see this, note first that 
if the measurement is unbiased the vacuum probability 
distribution Q{A, B) will be independent of the phase. 
Second, for the squeezed state |a,^), £,a* /a is indepen- 
dent of the phase arg(Q:). This in turn means that BA* /A 
is independent of the phase. Since 



Aa* = (1 + BA*/Ay 



l-\B[ 



(2.11) 



Aa* and therefore Q'{A, B) are independent of the phase. 

Since the probability distribution for the phase de- 
pends on the inner product between the two states, the 
variance in the measured phase will approximately be the 
sum of the intrinsic phase variance and the phase vari- 
ance of the squeezed state la,^). The maximum over- 
lap between the states will be when the squeezed state 
has about the same photon number as the input state. 
This means that the theoretical limit to the phase vari- 
ance that is introduced by the measurement is the phase 
variance of the squeezed state that has the same photon 
number as the input state and has been optimized for 
minimum intrinsic phase variance. Since the phase vari- 
ance of a squeezed state optimized for minimum intrinsic 
phase variance is \nn/ {An^) in the limit of large n 
this is also the limit to the introduced phase variance. 

The photon number of the squeezed state at maximum 
overlap will be mainly determined by the photon number 
of the input, but the degree and direction of squeezing 
(parametrized by ^) will be determined by the multi- 
plying factor Q'{A,B). The multiplying factor can be 
expressed as a function of n and for which we will 
use the same symbol Q' , even though it is a new func- 
tion Q'{n, C). Here n is the mean photon number for the 
state \a,£) (and will be close to the photon number 
of the input state), and C = ^a*/a is ^ with the phase of 
a scaled out. The multiplying factor will tend to be con- 
centrated along a particular line, effectively giving C, as 
a function of n. In order to obtain the theoretical limit, 
the measurement scheme must give a multiplying factor 
Q'(n, C) that tends to give values of C, for each n that are 
the same as for optimized squeezed states. 

We can determine the approximate variation of Q with 
n in the multiplying factor if we can estimate how it varies 
for measurements on a coherent state. If we consider 
measurements on a coherent state with real amplitude 
a"^, then the maximum overlap with the state |a,^) will 
be for a^ ~ a. We use a^ without a subscript to indicate 
the initial coherent amplitude before the measurement. 

If we are using an adaptive scheme with intermediate 
phase estimates that are unbiased, it is easy to see that 
the maximum probability will be for B real and therefore 
also A real. These results imply that 



a : 



A{l + B) 

1-B2 



A 



1 -B 



(2.12) 
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In turn this gives C as 

C w -atanh(l - Aja) 



^1 ^ 



2 2VI 



(2.13) 
(2.14) 



Since the value of C is governed by the multiplying factor 
Q'{n,C), this result for ( should hold for more general 
input states. 

From Rcf. ^ the phase variance of a squeezed state is 



"0 + 1 



(2.15) 



where no = ne^'^ for real C. This is minimized asymptot- 
ically as 



Inn + A 

4n2 ' 



where A « 2.43, for 



no « ln(4n) - - ln(27r). 



(2.16) 



(2.17) 



If we use the result obtained for C in Eq- ( 2.14 ) we find 
that 



no 



(2.18) 



This result means that in order for the measurement to 
be optimal, \A\ should scale with n as 



1^1 oc 



Inn 



(2.19) 



For the case of mark II measurements we have the result 
that |A| = 1 1^, which is why these measurements are 
not optimal. Note that if we substitute \A\ = 1 into 
the expression (2.18) to find no, and substitute that into 
Eq. ( p.l5| ), we obtain the correct result for the mark II 
introduced phase variance, 



(A^2 



-1.5 



(2.20) 



III. IMPROVED FEEDBACK 



Now we have the result that for optimal feedback |^| 
should decrease with photon number. Therefore in order 
to improve the phase measurement scheme we want one 
that gives |^| < 1. To see in general how this can be 
achieved, consider a coherent state with amplitude 
and determine the Ito SDE (stochastic differential equa- 
tion) for 

d\A,\^ = AlidA,) + {dA^A, + (dADidA,) (3.1) 
= A;e**(^)/(z;)dw + e-'*('')/(t;)dwA„ + dv (3.2) 
= [|^.|/(w)2Re(e'*('')e-*^^) + l]dv, (3.3) 



where if^ = aigAy. In terms of the phase estimate ify 
^{v) — Tr/2 this becomes 



d\A,\^ = [1 + 2\A,\I{v) sin((^;^ - ip,)]dv. 



(3.4) 



If we take the expectation value of I(v) and simplify we 
get 

(/(i;)) = -2|a^|sin(^, ~^), (3.5) 

where (p = arga'^. If we use this result the expectation 
value for the increment in \Ay\'^ is 

((i|A„p) = [1 - 4|A„||a'^| sin(<^,„ - ip) sin((^^^ - ip.^)] dv. 

(3.6) 

The first term on its own will give \A\ = 1, and in order 
to get < 1 the two sines must have the same sign. 
This will be the case if the phase estimate is between the 
actual phase and the phase of Ay. It is for this reason 
that we consider phase estimates that are intermediate 
between the best phase estimate and the phase of Ay, 
i.e., of the form 

^{v) = arg[(A.^; + ByA:y-<^^A<% (3.7) 

In general, smaller values of |A| can be obtained by using 
smaller values of e. This is because (py tends to be a worse 
phas e es timate, thus making it possible for the sines in 
Eq. (3.6) to be larger. Note that it is far too simplistic 
to use the best phase estimate (i.e., with e = 0), as we 
need to adjust e in order to make no closer to optimal. 



IV. SIMULATION METHOD 

The easiest input states to use for numerical simula- 
tions are coherent states, as they remain coherent with 
a deterministically decaying amplitude. However, in or- 
der to estimate the phase variance that is introduced by 
the measurement this would be very inefficient, as the 
phase variance would be dominated by the intrinsic phase 
variance. It is almost as easy (and much more efficient) 
to perform calculations on squeezed states, as squeezed 
states remain squeezed states under the stochastic evolu- 
tion, and only the two squeezing parameters need be kept 
track of. The best squeezed states to use are those op- 
timized for minimum intrinsic phase variance. For these 
states the total phase variance will be approximately 
twice the intrinsic phase variance when the measurements 
are close to optimal. 

To determine the SDE's for the squeezing parameters, 
we must first consider the SDE for the state. For dyne de- 
tection the stochastic evolution of the conditioned state 
vector is 

~ 2" 



d\m) = 



dt 



(a^j + 7*a) 



+dN{t) 



ae 



l7l 



+7*)(a + 7)> 
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where a is the annihilation operator for the mode, I7I ^ 1 
is the ampUtude of the local oscillator, and $ = arg 7 is 
its phase. Here the mode being measured is assumed to 
come from a cavity with an intensity decay rate equal to 
unity. The point process dN{t) has a mean ndt, where 



(4.2) 



The equation given in |^] differs from Eq. (4.1) by a trivial 
phase factor. The form above is given because it is not 
possible to directly take the limit of large local oscillator 
amplitude using the form given in Q. To take the limit 
of large local oscillator amplitude we approximate the 
Poisson process SN{t) by a Gaussian process 



SN{t) « KSt + y/^SW{t), 



(4.3) 



where dW{t) is a Gaussian random variable of zero mean 
and variance St. Then we find that in the limit of large 
I7I we have 



d^it)) = + axe"** - xV2)dt 



(4.4) 



where 



-i<I> 



(at)e**). (4.5) 



In order to determine the SDE's for the squeezing pa- 
rameters, we use the method of Rigo et al. [|lO|. Squeezed 
states obey the relation 



(4.6) 



The squeezing parameters and Bf are related to the 
usual squeezin g pa rameters i n th e same way as A and 
B arc in Eq. (2^) and Eq. (2/7). In the Stratonovich 
formalism 



(« - at - Af)dm)) 



{dB^a) 



■dA^,mt))- (4.7) 



Converting the SDE for the state to the Stratonovich 
form in the usual way we find 



dm)) 



+ (ae 
1 

~2 



2 

— 1$ 



2 

x)dW 
adie""^) - dx]dW 



2axe 



\m)- 



dt 



(4.8) 



Here we have included the increments d(e"**) and dx be- 
cause the phase of the local oscillator can vary stochasti- 
cally. Using this form of the equation, the left hand side 
of Eq. (4.7) evaluates to 



dt[-ia^B^+Af/2)~BfiB^a^+Af)e-'^^+2Bfxe 



+dW 



Bf e-* - i^f d(e-*)] I 



\m)- 



This gives us the SDE's for the squeezing parameters. 



dB: 



^Bf(l + e-2'*Bf)dt, 



(4.9) 



dAf = -^Afdt+B^{a^){l + e-^'^B^)dt+B^e-'^dW 



-B^d{e-^'')dW. 



(4.10) 



From these we find that the Stratonovich SDE for the 
standard (nonscaled) amplitude af is 



da? = - ^afdt 



1 

1 BfdW 



2 1 



[(Bn*d(e^*) + rf(e-*)]. 



(4.11) 



Gonverting back to the Ito SDE, we get 
1 



dat = -- 



'dt 



1 



(4.12) 



The SDE for Bf is unchanged under the change to Ito 
form. If we take the signal to be I{t)5t = {5N+ - 6N-)/l3 
(for consistency with Ref. [||), then take the limit of large 
oscillator amplitude and small time intervals St, we ob- 
tain 

I{t)dt = 2Re(af e"'*(*))d^ -I- dW{t). (4.13) 
The parameters At and Bt are then defined as in w] by 



At = / e**e-^/2/(s)ds, 
Jo 

Jo 



e^'^e-'ds. 



(4.14) 
(4.15) 



In order to get rid of the exponential factors, we change 
the time variable to 



1 



(4.16) 



and we redefine the amplitude to remove the systematic 
variation: 



(4.17) 



Here we use the v subscript to indicate the scaled ampli- 
tude, and the t subscript to indicate the original, unsealed 
amplitude. Since these are equal to each other at zero 
time, there is no ambiguity in the initial amplitude . 
Reverting to our original definition of the signal (1.1), we 
find 



I[v)dv = 2Re(a^e-**(''))dv + dW[v). 



(4.18) 
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With these changes of variables, the definitions for 
and B„ become 



B,, 



e^'^I[u)du, 



(4.19) 
(4.20) 



The differential equations for the squeezing parameters 
become 



daZ = 



dv 

l-v 



1 B^dW{v) 

l-v 1-|R^|2 



(4.21) 
(4.22) 



Initial calculations were performed using these equations, 
but there is a further simplification that can be made. 
The solution for B^. is 



B: 



l-v 



{Bl)-^-Bl 



(4.23) 



For calculations with time-dependent e this solution for 
B^ was used rather than solving a separate differential 
equation for B^ . 



V. RESULTS 

First we will describe the results for constant e. For 
each mean photon number, e was varied to find the value 
that gave the minimum phase variance. This method 
does not give results close to the theoretical limit for 
photon numbers above about 5000, but the phase vari- 
ances continue to get smaller as compared to the phase 
variances for mark II measurements. This indicates that 
the results are following a different scaling law, and fit- 
ting techniques give the power for the introduced phase 
variance as 1.685 ± 0.007. The data and the fitted line 
along with the heterodyne and mark II cases and the 
theoretical limit are shown in Fig ^ These results are 
a significant improvement over the mark II case, but are 
still significantly above the theoretical limit. 

In order to improve on this result we must vary e during 
the measurement. The value of e that we found to give 
the best result was 



e{v) 



\B„ 



1 



(5.1) 



The reason for the multiplying factor of (t;^ — |i3t,p)/|Ci,| 
is that it is an estimator for l/|a'^|. This means that the 
value of e tends to be smaller for larger photon num- 
bers, resulting in smaller values of \A\. The reason for 
the factor of ^w/(l — v) is that it makes the value of e 
close to zero initially, and very large near the end of the 
measurement. 

This second factor was found essentially by trial and er- 
ror, and is thought to be related to the fact that the phase 



Mo"' 
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FIG. 2: Phase variance for phase measurements with a con- 
stant value of e plotted as a function of the photon number 
of the input state. The crosses are the values obtained by 
stochastic integration and the continuous line is the fitted 
line. For comparison we have also plotted, in order from top 
to bottom, the variance for heterodyne measurements (dashed 
line), for mark II measurements (dash-dotted line), and the 
theoretical limit (dotted line). 

of varies stochastically during the measurement. Re- 
call that during the measurement we want the phase es- 
timate to be between the phase of and the phase of 
Ay. We only have an estimate of the phase of (the 
initial phase), so if we use a phase estimate that is too 
close to the actual phase when the phase variance of 
is large, the phase estimate is likely to be outside the 
interval between the phase of and the phase of Ay. 
Since the phase variance of increases with time, the 
value of e is increased as well, to prevent this happening. 
The results for this method are shown in Fig. as 



1.05 




10 10 
Photon number 

FIG. 3: Phase variance for phase measurements with a time- 
dependent e plotted as a function of the photon number of the 
input state. The phase variance is plotted as a ratio to the 
theoretical minimum phase variance (i.e., twice the intrinsic 
phase variance). 
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a ratio to the theoretical limit. As this shows, the re- 
sults are very close to the theoretical limit, and even for 
the largest photon number for which calculations have 
been performed the phase uncertainty is only about 4% 
above the theoretical limit. For these calculations the 
time steps used were approximately 



Av = 



25 



(5.2) 



where (A(/)2)th is the theoretical limit to the phase un- 
certainty. With these time steps the uncertainty due to 
the finite step size is approximately 1%. 

If the integration time step is reduced, while keeping 
the time interval at which the phase estimates are up- 
dated constant, the phase variance converges. If, how- 
ever, the phase estimates are updated at smaller and 
smaller time intervals then the phase variance does not 
converge. For example, the phase uncertainty for mea- 
surements on an optimized squeezed state with a photon 
number of 1577 is 1.54 x 10~^ if we use the time steps 
given above. If, however, we use time steps that are 100 
times smaller, then the phase variance is 1.93 x 10~^, 
and if the time steps are 1000 times smaller the phase 
variance is 2.13 x 10~^. These results indicate that the 
phase estimates must be incremented in finite time in- 
tervals for this method to give good results, and the size 
of the time steps that should be used depends on the 
photon number. The phase variance is not strongly de- 
pendent on these time steps, however, and only an order 
of magnitude estimate of the photon number is required. 



VI. EVALUATION OF METHOD 

A problem with determing the phase variance by the 
method above is that, for highly squeezed states (that are 
close to optimized for minimum phase variance), a signif- 
icant contribution to the phase variance is from low prob- 
ability results around tt. In obtaining numerical results 
the actual phase variance for the measurement will tend 
to be underestimated because the results from around tt 
are obtained too rarely for good statistics. It would re- 
quire an extremely large number of samples to estimate 
this contribution. However, we can estimate it nonsta- 
tistically as follows. 

Recall that in order to have a measurement that is 
close to optimum the multiplying factor Q'{n,C,) should 
give values of ( for each n that are close to optimized for 
minimum phase uncertainty. To test this for the phase 
measurement scheme described above, the n and ^ were 
determined from the values of A and B from the samples. 
The resulting data along with the line for optimized C are 
plotted in Fig. ^. The imaginary part of C should be zero 
for optimum measurements, and is small for these re- 
sults. Therefore in Fig. ^ we have plotted the real part 
(h- As can be seen, the vast majority of the data points 
are below the line, indicating greater squeezing than op- 
timum. This means that if the low probability results 
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FIG. 4: Values of <^_r and n (calculated from A and B) result- 
ing from measurements on squeezed states of various mean 
photon numbers. The variation of ^ with n for optimum 
squeezed states is also plotted (continuous line). 

around tt are taken into account the phase variance for 
these measurements will be above the theoretical limit. 

First we consider the effect of variations in the modulus 
of leaving consideration of error in the phase till later. 
In order to estimate how far above the theoretical limit 
the actual phase variance is, we make a quadratic approx- 
imation to the expression for the phase variance. From 
Q the expression for the phase variance of a squeezed 
state is, for real 



in 



^2erfc(V2^e^). 
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Taking the derivative with respect to C gives 




(6.1) 
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Taking the second derivative and using the fact that 
the expression above is zero for minimum phase variance 
gives 



^(1 



4no). 
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This means that for values of ^ close to optimum the 
increase in the phase variance over the optimum value is 



A( 



4nn 



(6.4) 



The main contribution to the phase uncertainty is 
no/(4n2), so the increase in the phase uncertainty as a 
ratio to the minimum phase uncertainty is 



(A|Cl)'(l + 4no) 



(6.5) 



This estimate indicates that the actual phase variance 
for the measurement scheme described above can be sig- 
nificantly larger than the intrinsic phase variance. For 
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example, for a mean photon number of about 332 000 
the rms 

deviation of |C| from the optimum value is only about 
0.16, but a squeezed state with |^| differing this much 
from optimum will have a phase variance more than twice 
the optimum value. This indicates that if the low prob- 
ability results around tt are taken into account the in- 
troduced phase variance is actually more than twice the 
theoretical limit. 

Next, we estimate the contribution from error in the 
phase (rather than the modulus) of C- For a squeezed 
state with real a the intrinsic uncertainty in the zero 
quadrature is 



cos — 
2 



2|C| • 2 M 

e sm — , 



where /i = arg^. Since = 2asin(0) 
trinsic uncertainty in the phase is 



2aq 



e-^l'^lcos^(Af/2)-Ke^l'^lsin^(^/2) 
4fi 



(6.6) 
the in- 

(6.7) 



If the phase of C, is small, we can make the approximation 



-2|CI+e^lCI^2^4 
4n 



(6.8) 



Clearly the first term in the numerator is just the original 
phase variance, and the second term is the excess phase 
variance due to the error in the phase of C,. Therefore 
the extra phase variance due to error in the phase of C, is 
given by 



(AargQ' 



(6.9) 



Using this estimate on the previous example it can 
be seen that this is not so much of a problem, with the 
introduced phase uncertainty being increased by less than 
3% by this factor. 



VII. IMPROVED METHOD 

The problem of the large contribution of the low prob- 
ability results around tt can be effectively eliminated in 
the following way. At each time step the photon num- 
ber is estimated from the values of and _B„, and the 
optimum value of Q is estimated using the asymptotic 
formula in Then if C,r (the real part of Q is too far 
below the optimum value, rather than using the feedback 
phase above, we use 



*(^) = \ arg 



By - w^tanhlCoptl 



(7.1) 



Using this feedback phase takes By directly towards the 
optimum value. To see this, note that the optimum value 
of By is 

Cy 



V— tanhlCoptI 



(7.2) 
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FIG. 5: Contributions to the phase uncertainty from error 
in the magnitude of ^ (continuous line) and the phase of 
(dash-dotted line). Dividing factors of 1.1 and 1.2 are used 
for photon numbers of 5 x lO'' and 2 x 10^, respectively. For 
these photon numbers A = 5 x 10~*, otherwise A — 10^'^. 
The contributions are plotted as a ratio to the theoretical 
minimum introduced phase uncertainty. 

Taki ng t he exponential of the feedback phase given by 
Eq. ^ gives e^** (x By - B°p\ so dBy cx B°p* - By. 

The details of exactly when is considered too far be- 
low optimum can be varied endlessly, but for the results 
that will be presented here we use this alternate phase 
estimate after time v = 0.9 and when 



ICI > |Copt|e^l""l'(^-''\ 



(7.3) 



where Copt is the estimated optimum value of C and ay 
is Cy/{v'^ — \By\'^). Using the exponential multiplying 
factor means that the alternative feedback is used only 
toward the end of the measurement. Only considering the 
alternative feedback in the last 10% of the measurement 
is necessary for the smaller photon numbers, where Eq. 
(7.3) is too weak a restriction. 

Another variation from the previous scheme is that, for 
the larger photon numbers, the values of e given by the 
original expression were reduced. The above correction 
corrects only for values of (^r that are below optimum, 
and for the larger photon numbers many of the uncor- 
rected values of tend to be above optimum (see Fig. 
^). The corrections will still work well, however, if we use 
a dividing factor to bring the uncorrected values below 
the line. For the second largest photon number tested of 
around 5 x 10^, the best re sult s were obtained when the 
values of e as given by Eq. (5.1) were divided by 1.1. For 
the largest mean photon number tested, 2 x 10^, the best 
results were obtained for a dividing factor of 1.2. The 
value of A that gave the best results with these dividing 
factors was 5 x 10"''. For all other mean photon numbers 
the value of A used was 10"'^. 

The estimated contributions to the phase variance due 
to error in the magnitude and phase of C, are plotted in 
Fig. |. As can be seen, the contribution due to error in 
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the magnitude of C is very small, around 1.5% for the 
larger photon numbers tested. The contribution due to 
the error in the phase of C is a bit larger, but it still does 
not rise above 3%. Thus we can see that the introduced 
phase variance can be made very close to the theoretical 
limit, within 5% for the largest photon number tested. 

With this modified technique the phase variance again 
does not converge as the feedback phase is updated in 
smaller and smaller time intervals. The phase variance is 
less dependent on the time step with this technique, how- 
ever. For example, for a mean photon number of 1577 
the total phase variance for measurements on an opti- 
mized squeezed state increases by only about 9% as the 
time steps are reduced by a factor of 1000. In contrast, 
the phase variance increases by a factor of 38% for the 
previous technique. 



VIII. CONCLUSIONS 

Any estimate of an initially unknown optical phase 
made using standard devices (linear optical and opto- 
electronic devices, a local oscillator, and photodetectors) 
must have an uncertainty above the intrinsic quantum 
uncertainty in the phase of the input state. The mini- 
mum magnitude of the added phase variance was deter- 
mined in Ref. H to scale asymptotically as 



Ini 



4(f 



(8.1) 



where is the mean photon number of the input state. 
Previous phase measurement schemes do not approach 
this theoretical limit. In this paper we have shown that 
an adaptive phase measurement scheme not previously 
considered can attain this theoretical limit. In other 
words, we have determined what is essentially the best 
possible phase measurement technique. 



In practice, phase measurements are currently limited 
by detector inefficiency. For detector efficiency 7] the in- 
troduced phase variance cannot be reduced below ||| 



■2) 



When the mark II phase variance is less than this there 
is not likely to be any significant advantage to using a 
more advanced feedback scheme. For the best photode- 
tectors available today, with around 98% efficiency [ pT| , 
the mark II phase variance falls below this limit for pho- 
ton numbers above 1000. Below this photon number the 
mark II phase variance is never more than about 27% 
above the limits determined using Eqs. (8.1) and (^.2|), 
so only relatively small improvements can be obtained by 
using a more advanced feedback scheme. 

Nevertheless, the technology is always improving, and 
there is no fundamental reason why photodetectors can- 
not be built with efficiencies extremely close to 1 [ p^ . 
When very efficient photodetectors are developed, the 
feedback techniques described here have the potential to 
give great improvements in the accuracy of phase mea- 
surements for applications where there is a limitation on 
the photon number that can be used. The other detri- 
mental factors are relatively minor, although the time de- 
lay in the feedback loop will become significant for very 
short pulses. 

The primary significance of the result obtained in this 
paper is theoretical, however, as it represents the culmi- 
nation of the search for the best optical phase measure- 
ment schemes using standard devices. To do any better 
would require using nonlinear optical devices. For exam- 
ple, it is conceivable that down-converting some portion 
of the signal field, and then measuring the phase of the 
down-converted light, could enable the above theoretical 
limit to be surpassed. This is a question for future work. 
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